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Abstract Quasi one dimensional Bose-Einstein con- 
densates (BECs) in elongated traps exhibit significant 
phase fluctuations even at very low temperatures. We 
present recent experimental results on the dynamic trans- 
formation of phase fluctuations into density modulations 
during time-of-flight and show the excellent quantitative 
agreement with the theoretical prediction. In addition 
we confirm that under our experimental conditions, in 
the magnetic trap density modulations are strongly sup- 
pressed even when the phase fluctuates. The paper also 
discusses our theoretical results on control of the conden- 
sate phase by employing a time-dependent perturbation. 
Our results set important limitations on future applica- 
tions of BEG in precision atom interferometry and atom 
optics, but at the same time suggest pathways to over- 
come these limitations. 

PACS: 03.75.Fi, 32.80.Pj, 05.30.Jp 



1 Introduction 

Since the first experimental observation of Bose-Einstein 
condensates |^ a considerable effort has been devoted to 
studies of low dimensional ultra-cold trapped gases Q . 
There are several fundamental reasons for the interest 
in these systems. Firstly, atom-atom interactions can be 
modified in low dimensions Secondly, the physics 
in low dimensional systems can be remarkably differ- 
ent from the situation in three dimensions. For example, 
a one dimensional condensate in the limit of ultra-low 
densities should behave as a strongly correlated system 
of impenetrable bosons, the so-called Tonks gas The 
third, and perhaps the most important reason to study 
low temperature ID condensates, is the particular po- 
tential that such condensates offer for applications. Pro- 
posals of atom optics and precision atom interferometry 
with coherent matter waves rely on placing the conden- 
sate in a ID waveguide. Such a waveguide can be formed 



in appropriately designed traps Q or on the surface of 
an atom chip [|[|^ . One could then use the BEG on such 
a chip to perform beam splitting and, finally, beam re- 
combination in a waveguide. 

However, it has been shown that BEG in quasi ID 
geometries exhibit significant phase fiuctuations l^],^. 
Even at very low temperatures (of the order of fractions 
of the critical temperature Tc) the coherence length of 
such a phase-fluctuating BEG may be shorter than the 
size of the condensate. The condensate then exhibits a 
spatially varying phase pattern and is commonly called 
a quasicondensate. In the Thomas- Fermi regime [^o|, 
that is, when the nonlinear mean field energy is much 
larger than the kinetic energy, density fiuctuations are 
suppressed due to their energetic cost. That is not the 
case, however, for the phase fiuctuations. The prediction 
of phase fluctuations |^ was experimentally confirmed 
in |ll|,0. 

The results of our experiments set important limi- 
tations on the future applications of BEG in interfer- 
ometry which employ ID waveguides. It is, therefore, 
particularly important to study and characterize phase 
fluctuations in elongated BECs and, even more, to de- 
velop methods that allow to control the phase of these 
condensate. 

The first part of this paper concerns the experimental 
investigation of two aspects of phase-fluctuating BECs. 
First the dynamical behavior during a time-of-flight after 
release from a strongly elongated trap is examined. Very 
good quantitative agreement between theory and exper- 
iment is established for these measurements. A second 
set of measurements is devoted to the investigation of 
the density of a phase-fluctuating trapped condensate. 
We conflrm that the density fluctuations are strongly 
suppressed and show that the second order correlation 
function ^'^■'(O) is largely independent of the amount of 
phase fluctuation in our experimental regime. 

The second part of the paper deals with theoretical 
approaches to overcome the effects of phase fluctuations. 
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We propose to use the technique of parametric resonance 
to control the overaU phase of the condensate. It should 
be noted that parametric resonance has already been 
discussed in the context of BEC but those 

studies were aimed at employing the parametric reso- 
nance to create patterns and textures in the density of 
trapped BECs. Clearly, for the reasons mentioned above, 
such a goal is not easily accomplished in highly nonlinear 
systems for which fluctuations of the mean field energy 
are costly. We, however, aim at controlling the phase of a 
quasi ID condensate by applying the standard technique 
of modulating the trap frequency [p^ . Thus, rather than 
to excite density modulations of the condensate, we in- 
tend to engineer its phase, such that some modes con- 
tributing to the phase fluctuations can be suppressed or 
enhanced without destroying the condensate. 

The paper is organized as follows: In Section 2 we 
briefly review the physics of phase fluctuations in elon- 
gated traps, and recall the expressions describing the 
transformation of phase fluctuations into density modu- 
lations [pl| . In Section 3 we present recent experimental 
results which expand the observations reported in our 
earlier work. These results show that a full quantitative 
characterization of phase fluctuations can be achieved 
and that density modulations are suppressed even when 
phase fluctuations are present in a trapped condensate. 
In Section 4, we describe the method of employing para- 
metric resonance to enhance or suppress selected modes 
of the phase fluctuations. We present model calculations 
based on numerical solutions of the ID Gross-Pitaevskii 
(GP) equation which governs the dynamics of the con- 
densate. 



2 Description of the phase fluctuations 

For BECs in 3D trapping geometries the fluctuations of 
density and phase are only important in a narrow tem- 
perature range near the BEC transition temperature Tq 
(critical fluctuations). Outside this region, the fluctua- 
tions are suppressed and the condensate is phase and 
density coherent. In reduced dimensions, the situation 
can be rather different (see and refs. therein): the 

axial phase fluctuations can manifest themselves even at 
temperatures far below Tq. 

In the following we consider a cylindrically symmet- 
ric condensate in the Thomas-Fermi regime, where the 
mean field (repulsive) interparticle interaction greatly 
exceeds the radial {hcup) and axial (hLOx) trap energies. 
At T = its density profile has the well-known shape 
no (p, x) = nom ( 1 - / RIf - /^tf ) > where nom = m/s 
is the maximum condensate density, fi is the chemical 
potential, g = Anfi^a/m is the interaction coupling con- 
stant, m is the atomic mass, and a > the scattering 
length. Under the condition LUp ^ lu^, the radial size 
of the condensate, given by the Thomas-Fermi radius 
i?TF = (2Ai/TOWp)^/^, is much smaller than the axial 



size, given by the corresponding Thomas-Fermi length 
L^p = (2M/mL^2)i/2_ 

The phase fluctuations can be described by solving 
the Bogoliubov-de Gcnnes equations describing elemen- 
tary excitations of the condensate. The quantum field 
annihilation operator of atoms can be written as V'(r) = 
y/nQ{r) exp(j(/)(r)), where 0(r) is the operator of the 
phase, and the density fluctuations have already been 
neglected following the arguments discussed above. The 
single-particle correlation function is then expressed 
through the mean square fluctuations of the phase (see 
e.g. iisl): 



(^t(r)V-(r')) -v/«o(r)no(r')exp{-([<50(r, r')]2)/2}, 

with S(j){r, r') = 0(r) — (pir'). The operator 0(r) is given 
by (see e.g. [|9|) 



(r) = [4no{r] 



]-V2^/+(r)a^.+h.c., 



(2) 



where dj is the annihilation operator of the excitation 
with quantum number j and energy Cj , /^^ — Uj + Vj , 
and the excitation functions Uj , Vj are determined by the 
Bogoliubov-de Gennes equations in the Thomas-Fermi 
limit. 

The "low-energy" axial excitations (with energies 
ej < hujp) have wavelengths larger than i?TF and ex- 
hibit a pronounced ID behavior. Hence, one expects 
that these excitations give the most important contri- 
bution to the long-wave axial fluctuations of the phase. 
The solution of the hydrodynamic counterpart of the 
Bogoliubov-de Gennes equations in the Thomas-Fermi 
limit for such low-energy axial modes yields the spec- 
trum €j — fiLOx yj j {j + 3) /4 j2^, where j is a positive 
integer. The wavefunctions f^' of these modes have the 
form 



/+(r) 



' (j + 2)(2j+3)gno(r) (1,1) 



(3) 



where -Pj'^'^'' are Jacobi polynomials. Note that the spe- 
cific form of the excitation spectrum and the x-depen- 
dence of the mode functions results in the quasi ID 
limit from the integration of the mode functions over 
the transverse directions, as pointed out for instance in 
Ref. @]. 

The thermal fiuctuations of the phase behave for dis- 
tances shorter than ~ 0.4 



{mx,x')r)T^si\x~xyL, 



(4) 



where is a measure for the amount of phase fluctua- 
tions present in the condensate and given by H 



Sl{T) = T/T^ = 32fikBT/15Noihu 



(5) 



where fee is the Boltzmann constant and = L^^-p/Sj^ 
can be interpreted as a phase coherence length. Sub- 
stituting the chemical potential in the Thomas-Fermi 
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regime in Eq. one obtains the direct dependence on 
the experimental parameters: 



Ltf/I^, = 16 



Vl53/2?i3 



2/5 



AT, 



3/5 



4/5 



(6) 



Strong phase fluctuations are thus associated with high 
trap aspect ratios A = LOp/iOxi weak axial confinement 
ujx, high temperatures T and small numbers of con- 
densed atom A^o- 

The signature of a fluctuating phase can be observed 
experimentally as density modulations (stripes) in the 
ballistic expansion. The formation of these stripes can 
be understood as follows: After switching off the trap, 
the mean field interaction rapidly decreases and the axial 
phase pattern is converted into a velocity distribution. 
This spatial velocity distribution leads to the appear- 
ance of stripes in the density after a time-of-flight. These 
stripes have been described analytically ||l2| using a gen- 
eralization of the Bogoliubov-de Gennes equation for the 
self-similar solution (see e.g. jl^,^,^). This equation 
describes the ballistic expansion in presence of fluctua- 
tions. By denoting the density as ng -I- 5n and the phase 
as 00 + we obtain the analytic expression for the rela- 
tive density fluctuations 



5fi 
no 



2E- 



-{<ij/hu}p) 



sm 



Ai?iC.p(l - (^)2) 



0.-, (7) 



where the sum extends over the axial modes, r — tOpt, 
t is the time-of-flight and (j)j is the contribution of the 
j-th mode to the phase operator in Eq. (||). 

By taking the thermal average of the square of Eq. 
(^ , one obtains a closed relation for the mean square 
density fluctuations psl 



\ = |^C^(^o,-p,o..,x,t), (8) 



where 



C^(iVo, Wp, w^,a;,t) = 

2( 0- + 3/2)^ 
2^- \Aa{l~{x/L^,Y) 



(9) 



with a = ^/huj^t. 

3 Experimental Results 

3.1 Experimental setup 



The experimental setup used to observe phase-fluctua- 
ting Bose-Einstein condensates has been described else- 
where Here we emphasize the main features of 



the apparatus and the technical changes that have led 
to the observation of more pronounced effects of phase 
fluctuations in the ballistic expansion of an elongated 
condensate. 

We have adapted the trap geometry to obtain aspect 
ratios A > 100 and reduced axial confinement. These 
high trap aspect ratios are achieved by axial decompres- 
sion of our cloverleaf trap p^ , lowering the currents in 
both the pinch coils and their bias coils. 

For the experiments reported here, a Bose-Einstein 
condensate of ®^Rb atoms in the \F =l,mp = —l) state 
is produced in the axially decompressed trap with axial 
and radial trapping frequencies 3.4 Hz and 380 Hz re- 
spectively, by using the following procedure. We load a 
magneto-optical trap with a few times 10^ atoms from 
a chirp slowed thermal beam, followed by a short period 
of sub-Doppler cooling. Part of the sub-Doppler cool- 
ing stage operates without repumping the atoms from 
the \F = 1) state, letting them accumulate in the lower 
hyperfine state. After optical pumping to the low field 
seeking mp — state, the atomic cloud is loaded into 
a cloverleaf magnetic trap. A cleaning light pulse on res- 
onance with the \F — 2) \F' = 2) transition then 
removes all remaining \F — 2) atoms from the trap. Fi- 
nally, the atomic cloud is adiabatically compressed to 
allow efficient radio frequency evaporative cooling. The 
condensate is loaded in the axially decompressed trap 
by using a two step evaporation procedure. The initial 
evaporative cooling is performed in our standard trap 
with aspect ratio 26, then the trap is ramped to the high 
aspect ratio configuration, and the final evaporation is 
carried out in that trapping geometry. After waiting for 
4 seconds to allow the system to reach an equilibrium 
state we switch off the trapping potential within 200 /iS 
and wait for a variable time-of-flight before detecting the 
atomic cloud by resonant absorption imaging. 

The experimental observation of phase fluctuations 
profits from this adapted trap configuration as indicated 
by Eq. (^) . Our present measurements show strong fluc- 
tuations both due to the high aspect ratio and the weak 
axial confinement of our trap. These lead to a large ax- 
ial condensate size (in the elongated trap, 2L ~ 300 /im) 
with respect to the detection resolution of ^ 5 /im, thus 
favoring the observation of structures in the condensate 
density after time-of-flight. 



3.2 Expansion dynamics 

Phase fluctuations transform into density modulations 
during ballistic expansion. It was shown that the gen- 
eral dependence on temperature, trapping potential and 
number of atoms given in Eq. (3) of [y agrees quali- 
tatively with the experimental observations. Most data 
was acquired after a long time-of-flight of 25 ms since 
the effect of phase fluctuations is best observed then. 
However, here we present observations of the dynamical 



4 



H. Kreutzmann et al. 



0,07 




4.6 ms 
9.6 ms 

14.6 ms 



19.6 ms 




- 200 200 
X [pm] 



Fig. 1 (a) Absorption images of the expanded clouds after 
various times-of- flight. The images show ~ 7 x 10* condensed 
atoms at a temperature of ~ 0.6 Tfc- The atomic clouds are 
released from a magnetic trap with Eixial and radial frequen- 
cies of 2tt X 3.4 Hz and 2ti x 380 Hz respectively. Due to the 
destructive imaging technique, each image shows a different 
condensate, (b) The line density profile of the leist time-of- 
flight image is compared with the corresponding bimodal fit- 
ting function. 



transformation of phase fluctuations into density modu- 
lations and quantitatively compare them to the theoret- 
ical predictions. 

Fig. |l|(a) shows typical images of the ballistically ex- 
panded clouds for various times-of-flight t. The atomic 
samples were released from a magnetic trap with axial 
and radial trapping frequencies of 3.4 Hz and 380 Hz, 
respectively. They contain ^ 7 x 10^ condensed atoms 
at a temperature of ~ 0.6 Tq. Under these experimental 
conditions, the phase coherence length in units of the 
condensate size is l,p/2L ^ 1/12. The usual anisotropic 
expansion of the condensates is clearly visible. The line 
density profile in Fig. |l|(b) shows the large influence of 



time-of-flight [ms] 

Fig. 2 Black circles: Measurements of the normalized spatial 
average of the phase fluctuations o"|ec ^ function of the 
time-of-flight. Each data point is the average of 15 measure- 
ments. Open squares: Theoretically predicted values. Line: 
Fit to the measured data as a guide to the eye. The error 
bars indicate the statistical error. 

phase fluctuations on the expansion of the condensate. 
To avoid high optical densities, which would influence 
the visibility of density modulations, we detuned the de- 
tection laser with respect to the atomic resonance for 
short times-of-flight. 

The amount of phase fluctuations present in a given 
BEC was determined by comparing the observed line 
density profile with the expected bimodal profile when 
no phase fluctuations are present. In this bimodal profile 
the BEC fraction is described by a parabolic Thomas- 
Fermi component and the thermal cloud by a Gaussian 
distribution, both integrated along the radial direction. 
The deviation of the experimental density profile from 
the fitted function is calculated at every data point in 
the interval —0.5 < x/L < 0.5 and normalized by the 
fitted condensate density at that point. The average of 
the squares of these values is called ct^ec ^"^^ represents a 
normalized spatial average of the phase fluctuations in a 
single condensate. The temperature and particle number 
of each condensate were determined by separate 2D fits 
to the absorption images. 

Note that phase fluctuations in an elongated BEC 
are stochastic. During the expansion, the instantaneous 
phase of the BEC at the time of release is converted into 
density modulations and therefore images taken under 
the same initial conditions can look significantly differ- 
ent. Indeed, we observe a large spread of our experi- 
mental data and, therefore, each data point in Fig. |^ 
represents the average of 15 measurements. 

Fig. H shows the observed density modulations (black 
circles) as a function of the free expansion time. All im- 
ages were taken under the conditions given in Fig. 0. The 
data clearly shows how the effect of phase fluctuations 
on the expansion dynamics increases with the time-of- 
flight. 
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To compare our measurements with theory, each ex- 
perimental realization (each condensate) was modeled 
individually as explained below and for each time-of- 
flight the average of all experimental realizations is com- 
pared to the average of all modeled realizations. Using 
the experimentally determined temperature and atom 
number, the expectation value of the density modula- 
tions for each realization was calculated according to 
Eq. (^). Analogous to the experimental procedure, a nor- 
malized spatial average over the range — 0.5<x/L<0.5 
was deduced. The theoretical prediction in Fig. || takes 
the limited resolution of our imaging system into ac- 
count m. 

We obtain very good quantitative agreement between 
our results and the theoretical predictions, confirming 
that the transfer mechanism of phase fluctuations into 
density modulations during time-of-flight is well under- 
stood. 



3. 3 Suppression of density fluctuations 

The theoretical description of phase fluctuations is based 
on the prediction that density modulations arc strong- 
ly suppressed in the magnetic trap. In the Thomas-Fermi 
regime, where the mean fleld energy dominates over the 
kinetic energy, the excitation of density modulations re- 
quires a high energetic cost on the order of the chemical 
potential. The measurements reported in Fig. |^ clearly 
show that the observed density modulations are sup- 
pressed as the time-of-flight approaches zero. However, 
this does not rule out the existence of density modula- 
tions on a length scale smaller than our experimental 
resolution. 

The suppression of density fluctuations can be veri- 
fied experimentally by measuring the second order cor- 
relation function ^'^^^(ri — of the fleld operator <f'(r). 
In second quantization formalism. 



(ri - ra) = 



n{vi)n{v2) 



(10) 



In particular, g'^'^^ (0) gives the correlation function of the 
atomic density for ri = r2. As shown in |Q , g^^^ (0) is di- 
rectly related to the expectation value of the interaction 
energy U by 



{U) 



(11) 



In the Thomas-Fermi regime, only the interaction energy 
(U) contributes to the kinetic energy of the condensate 
after ballistic expansion. The presence of density fluctu- 
ations would however lead to an increase of the observed 
release energy due to the presence of repulsive particle 
interactions. Hence, the ratio of the release energy and 
the calculated interaction energy in Thomas-Fermi ap- 
proximation gives the value of g^^^ (0). The energy due to 



2,0 
1,5 

s 1,0 

"bo 

0,5 
0,0 



—I ■ I ■ I ■ l_ 



100 200 300 400 500 
T[nK| 



2,5 
2,0 
1,5 
1,0 
0,5 
0,0 



Fig. 3 Open squares: Measurements of g'^'(O) as a function 
of temperature for phase-fluctuating \F— 1) condensates in a 
A = 30 trap. Black dots: Calculated phase coherence length 
Ig,, in units of the axial condensate size 2L, for the experi- 
mental conditions of each data point. The error bars indicate 
the maximum error. 



phase fluctuations is small compared to the interaction 
energy (for typical parameters ~ 0.5%) and can there- 
fore be neglected in release energy measurements. 

From earlier release-energy measurements, the prod- 
uct ag''^\0) has been determined |^,^|2j]. We have 
performed such a measurement for phase-fluctuating con- 
densates. 

Fig. |3| shows ^^^•'(O) (open squares) as a function of 
temperature for atoms in the \F = 1) state in a trap 
with A = 30. We use this weakly elongated trap, be- 
cause this trapping geometry enables us to use the sam- 
ple temperature to tune from nearly pure to strongly 
phase-fluctuating condensates experimentally. 

The same graph also shows the calculated phase co- 
herence length according to Equation (^ for the ex- 
perimental conditions of each data point (black circles). 
As the temperature of the sample is raised, the phase co- 
herence length decreases and becomes significantly shor- 
ter than the condensate size. Thus, our measurements 
clearly reach the regime of quasicondensation. However, 
the graph shows that the second-order correlation func- 
tion is largely independent of the condensate temper- 
ature T, and that it is consistent with the expected 
value g^^\0) = 1. These measurements give an upper 
limit for the density modulations present in a phase- 
fluctuating condensate, and clearly distinguish our sam- 
ples from thermal clouds where, due to bunching effects, 
g'^^^O) = 2 and large density fluctuations are present on 
small scales. 



4 Controlling the phase 

Parametric resonance refers to the exponentially large 
response of a system of a periodic external perturba- 
tion [M for some specific set of parameters. Although 
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parametric resonance is a very well established pheno- 
menon in linear systems, it was shown in psf that para- 
metric resonance can also occur in Bose-Einstein conden- 
sates. Normally one aims at using parametric resonance 
in order to modulate the condensate density by means of 
a relatively small periodic perturbation. However, in the 
Thomas-Fermi regime, as mentioned above, modulations 
of the density require large energies. In this Section we 
suggest a way to shape the overall phase by modulating 
the trap with a small amplitude for a short time. 



4-.1 Theoretical methods 

To illustrate the effect of the perturbation on the phase 
of a quasicondensate, we numerically solve the ID Gross- 
Pitaevskii equation for different initial conditions. We 
assume a zero temperature BEC in the Thomas-Fermi 
regime containing 5 x 10'' *^Rb atoms trapped in a ID 
confining potential with an (axial) frequency of lu^ = 
2tt X 14 Hz. The parabolic density profile is then given 
by no{x) — riom(l — I L'^w) where Ltf — (2/x/mw^)^/^ 
and TT-om = f'-l9- The coupling constant g in one di- 
mension can be derived by averaging the 3D interac- 
tions over the radial density profile. We first evolve the 
wavefunction according to the time dependent GP equa- 
tion in imaginary time to obtain the condensate wave- 
function at T = 0. We ensure that the initial phase 
of the condensate is constant before imposing a fluc- 
tuating phase corresponding to a fixed temperature T 
on the condensate wavefunction. In strictly ID, one ex- 
pects that for temperatures T^'^ T phase fluctu- 
ations are present whereas density fluctuations are sup- 
pressed |q|. H ere k^Td ~ Nhujx is the degeneracy tem- 
perature |]2^ , and Tif, = T^hcUx/lJ' is the characteristic 
temperature for the appearance of the phase fluctua- 
tions. 

Applying the same procedures for the ID case as 
described in Section 2 for the 3D case, the field an- 
nihilation operator can be written as iplx^t = 0) = 
y^rio{x) exu(i(l){x)). The operator of the phase is given 
by (see gg) 



1 



v4no(x) 



^/+(a:)a,+h.c. 



(12) 



where, as before, dj is the annihilation operator of the 
excitation with quantum number j, — Uj+Vj and uj, 
Vj are the excitation functions determined through the 
Bogoliubov-de Genncs equations in the Thomas-Fermi 
limit. The solution of these ID equations gives the spec- 
trum ej = hujx\/ j{j + l)/2 The functions have 
now the form 



where Pj (z) are the Legendre polynomials and z = x/Ltf- 
Note, that these solutions correspond to the strictly ID 
case and differ from those obtained in the quasi ID case 
described by (^). The qualitative and, to great extent, 
quantitative character of the solutions (^3|) and (^ is, 
however, similar. 

The random phase of the condensate is numerically 
simulated by replacing the operators (ij and in Eq. ( p^ ) 
by Gaussian random variables aj and a*, with the cor- 
relation {a*aj') — SjjiNj, where Nj is the occupation 
number for the quasiparticle mode j for a given temper- 
ature T. 

Once the phase is imposed, we apply a periodic per- 
turbation a sin(cL'si) to the trapping potential: 

uj^ ^li^ + Q! sin(wsi)). 



(14) 



'(j + 1/2) 



^(1 



Z2) 



(13) 



The amplitude a of the perturbation is small and ranges 
between 0.05 and 0.20 ^ 1, whereas the frequency ojs 
ranges from uj^ to several times Wx- The perturbation 
acts for a short period of time t (the simulations were 
performed for 2, 6 and 12 periods of modulation 2tt/ujs). 
After switching off the perturbation the amplitudes of 
the different modes participating in the fluctuations of 
the phase are calculated. The thermal fluctuations of 
the phase for the temperature range considered here are 
mostly provided by the low excitation modes (j < 50). 
Wc calculate the amplitude of these modes as a function 
of the frequency of the perturbation. 

4-2 Numerical results 

Our calculations show that the positive and negative 
frequency modes couple resonantly and that by appro- 
priately tuning the frequency of the perturbation one 
can selectively suppress some modes and enhance others. 
The enhancement, however, is not particularly strong, 
since we do not amplify any mode (except the second 
mode which arises directly due to the parametric modu- 
lation) much above its initial value provided by Eq. (p^. 
Note that the odd modes can only be excited if the sym- 
metry is broken, e.g. by the existence of phase fluctua- 
tions. 

Figure ^ displays the amplitudes of the first two modes 
for different initial temperatures as a function of the fre- 
quency ujs- The second mode, which corresponds to the 
density breathing mode, is strongly enhanced due to the 
perturbation of the trap. This mode amplitude can be 
calculated in the Thomas-Fermi regime by using a self 
similar solution of the GP equation HJ,^ describing 
the dynamics of the bare condensate that implies ap- 
pearance of a phase quadratic in x. Indeed, only for suf- 
ficiently large temperatures (T ~ 0.5 Tc); its frequency 
dependence starts to deviate significantly from its be- 
havior at r = 0. 

To show the effect of the fluctuating phase on the 
final amplitudes, we subtract the phase a pure conden- 
sate would have acquired under the same perturbation 
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Fig. 4 (a) Amplitude of the flrst and (b) second mode as a 
function of the modulation frequency ug. The perturbation 
was applied for t = 2 x 2n/u!s with an amplitude of a = 0.05. 
The solid line indicates the final amplitude starting from a 
pure condensate (T = 0); crosses correspond to a quasicon- 
densate at initial temperature of T = 0.01 Tc, open squares 
to T = 0.1 Tc and circles to T = 0.4 Tc. The same set of 
random coefficients were used in all the cases. 

from our results. Figure |^ displays the modulus square 
of such amplitudes normalized to their initial values for 
the first four modes. The amplitudes of the high excited 
modes (app. j > 4) are practically equal and neglibigly 
small compared to the lower modes therefore the lower 
modes dominate the dynamics of the phase. Only those 
modes show clear maxima and minima as a function of 
the modulation frequency luq. Resonances between dif- 
ferent modes may appear when 



in, ± 



(15) 



where 1?^ = ei/Ti. In order to observe resonances it is 
necessary that the energy differences between the modes 
correspond to multiples of the modulation frequency, 
otherwise the system behaves irregularly. In the ID case, 
the excitation modes are quite regularly separated by 
approximately oJx/V^ (~ 0.71 Wj^). However, in nonlin- 
ear systems, the resonance frequencies undergo shifts 
which can already be quite significant at small pertur- 
bations 13 . For the first mode (Fig. |), we find clear 
resonances at los/wx — 1-4, 2.1 and 4.0. One is tempted 
to attribute the first two resonances to the coupling with 
j = 5 and J = 4 mode. The second mode, which domi- 
nates the dynamics — since this mode is enhanced by the 
perturbation — shows three maxima at wg/w^ ~ 1.4, 2.1 
and 3.6. The analysis of maxima for the third and higher 
modes becomes quite complex. For the third mode, we 
find maxima at ujs/i^x — 1-4, 1.7, 2.4, 3.4 and 5.0. Start- 
ing from the fourth mode it becomes difRcult to associate 
the maxima with specific resonances. 

It is important to stress here that for the lower modes, 
the position of the maxima does not depend on the ran- 
dom generated phase. This fact clearly confirms that 
we are dealing here with resonances between different 
modes. By increasing the perturbation time the spec- 
trum changes, as expected, since this is a time dependent 
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Fig. 5 Modulus square of the exciting mode amplitudes nor- 
malized to their initial value at f = 0. The contribution cor- 
responding to T = has been previously subtracted, (a) (b) 
(c) and (d) correspond to amplitudes of the first second, third 
and fourth mode for the same parameters as Fig. W. 



Fig. 6 The dotted line shows a typical phase pattern for 
Lus ~ ^TT X 14 Hz, A = 51 and initial temperature T = 0.4 Tc, 
the solid line shows the final phase after applying the time 
dependent perturbation for a duration of two cycles with a = 
0.05 and ujs ~ 5.7 Ux- 



effect, and more peaks appear. However with the excep- 
tion of the second one, no enhancement of the modes is 
observed. By increasing the amplitude a of the pertur- 
bation, the second mode becomes much more enhanced, 
but again the amplitudes of the other modes (normalized 
to their initial values at t = 0) do not. In other words, 
we do not find an exponential parametric amplification 
by tuning the frequency to 2]7i/n where n is a positive 
integer. Nevertheless, the fact that the amplitudes of the 
lower modes display well defined maxima and minima al- 
lows to shape the overall phase of the quasicondensate. 
As an example, in Fig. ^ we show the initial fluctuat- 
ing phase of a quasicondensate at T = 0.4 Tc, and its 
final phase after the perturbation is switched off. In the 
range Ix/LtfI < 0.5 the final phase becomes practically 
constant demonstrating that the initial fluctuations can 
be significantly suppressed using the parametric modu- 
lation technique. To achieve this specific phase we chose 
the frequency tos such that the lower modes have ap- 
proximately the same amplitude value. 
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5 Conclusion and outlook 

In this paper we analyze fluctuations of the phase of 
quasi ID condensates at temperatures of the order of 
fractions of Tq- The dynamical transformation of phase 
fluctuations into density modulations was observed as 
function of free expansion time. A detailed comparison 
of the statistical average of the modulation with the 
theoretical prediction shows excellent agreement. It was 
also verified experimentally that density modulations of 
a phase fluctuating BEC in the trap are strongly sup- 
pressed. The quantitative understanding of the transfer 
of phase fluctuations into density modulations opens a 
pathway to use phase fluctuations for condensate ther- 
mometry. 

Phase fluctuations impose restrictions on the applica- 
bility of quasi ID condensates in atom optics and preci- 
sion atom interferometry. We show that a possible mech- 
anism to control the overall phase of a quasi ID conden- 
sate is to use a small periodic perturbation of the trap. In 
this way, only few of the modes responsible for the phase 
fluctuations become relevant. Furthermore, by properly 
tuning the frequency, time and amplitude of the exter- 
nal perturbation one can resonantly couple the different 
relevant low order modes, so that some of them are en- 
hanced and others are inhibited. This method suggest 
a pathway to control phase fluctuations and overcome 
their undesired effects. 

We acknowledge valuable discussions with G. Shlyap- 
nikov, D. S. Petrov and C. Borde. This work is sup- 
ported by the Deutsche Forschungsgemeinschaft within 
the SFB 407 and the Schwerpunktprogramm " Wechsel- 
wirkungen ultrakalter atomarer und molekularer Gase" , 
Alexander von Humboldt Stiftung, and the European 
Science Foundation (ESF) within the BEC2000-K pro- 
gramme. 
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